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In the inert scalar doublet framework, we analyze what would be the effect of a B — L asymmetry that could 
have been produced in the Universe thermal bath at high temperature. We show that, unless the “A.g” scalar 
interaction is tiny, this asymmetry is automatically reprocessed in part into an inert scalar asymmetry that could 
be at the origin of dark matter today. Along this scenario, the inert mass scale lies in the few-TeV range and 
direct detection constraints require that the inert scalar particles decay into a lighter dark matter particle which, 
as the inert doublet, is odd under a Z 2 symmetry. 


I. INTRODUCTION 


The similarity of baryonic and dark matter (DM) abun¬ 
dances, determined by the observation of the Cosmic Mi¬ 
crowave Background (CMB) anisotropies, CIdm/^b = 5.4± 
0.1 m, has motivated a long series of scenarios were both 
abundances have a related or even very same origin. Since 
the baryon asymmetry is to a very good approximation to¬ 
tally asymmetric - no primordial population of antibaryons 
has been observed in the Universe - a common origin of both 
abundances suggests that the DM abundance today would 
be associated to the generation of a DM particle-antiparticle 
asymmetry (see e.g. the reviews of Refs. EHSi). In the follow¬ 
ing, we will show how this can be realized from the generation 
of a scalar inert doublet asymmetry. The inert scalar doublet 
DM framework (IDM) fhHll simply consists in adding to the 
Standard Model (SM) a single scalar doublet, H 2 , odd under a 
Z 2 symmetry. The most general scalar potential is in this case 


V =m\\Hi\^+ml\H2\^+Xi\Hi\^ + X2\H2\‘^ 


+ X3\Hi\^\H2\^ + X4\hIH2\^ 


2 



+ h.c. 


( 1 ) 


where the SM and the inert scalar doublets can be written as 






and H 2 = 


( 2 ) 


with = {h + i(^^)/y/2 and r]*’ = {H^ + iA^)/y/2. In the 
scalar potential, is assumed to be positive to insure that 
H 2 doesn’t acquire a vev, so that it’s lightest (neutral) compo¬ 
nent is stable, unless there exists a lighter Z 2 odd particle into 
which it can decay. This is the possibility we will ultimately 
consider in the following, in order to satisfy the direct detec¬ 


tion constraints, see Section IV But before, up to the end of 
Section III, let’s restrict ourselves to the usual IDM setup and 
see what happens in this minimal framework. 

Prior to electroweak symmetry breaking (EWSB), all H 2 
components have mass mH 2 — m 2 , whereas after EWSB (v = 
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246 GeV), they get split in mass 

= 1712 Xf^ov^ , ni^o = 1712 -f A.^0 , iri^+ — tn^ -t- Xf^cv^ , 

(3) 

with Xf{c =X2,(2 and Xj^o = (Xg -|-A,4±A,5)/2. In the fol¬ 
lowing, we will assume, without loss of generality, that X 5 is 
negative, so that is lighter than A^. Various well-known 
constraints hold on the parameters of the theory. Tree level 
vacuum stability requires X\p > 0, Xfjo^^ojjc > —VX 1 X 2 « 
—Q36s/Xp,. EW precision test observables require AT ~ 
(m^+ —m^o) {m^+ —mfjo) /127t^av^ < 10^'. Z decay width 
constraint at LEP requires m^o -f > mz and > mz/2. 
Direct detection constraint importantly requires that the Z ex¬ 
change diagram is kinematically forbidden, i.e. — nifjo > 
A'rpDM/2^ where Pdmc is the DM halo velocity with respect 
to the earth, and pr = m^orng^j^/[lUfjo -\-mgf) is the reduced 
mass of the system for the nucleus used by the experi¬ 
ment. For nifjo ^ and Xenon nucleus, using an average 
velocity of ^ 270 km/s, this constraint can be rephrased as 
m^o — mjjo > ?>m,„i„ ^ 50 keV. Taking into account the veloc¬ 
ity distribution around this central value, and the recoil en¬ 
ergy sensitivity of the experiments, the minimum splitting be¬ 
comes 5mmm "^180 keV, although a more robust constraint is 
~ 100 keV HO), which translates as 

|X,|>3.3.1«-‘.(|a).(^). ,4, 

It is well known that the IDM can account for the observed 
DM relic abundance via the usual freeze-out mechanism, and 
be in agreement with direct detection constraints, for DM 
masses in the ranges ^ [50,80] GeV and above ^ 540 GeV, 
up to the ^ 40-50 TeV unitarity bound @ |7] [TTl [TJ]- As 
we will show, it could also be responsible for the DM relic 
density in an asymmetric way. 


II. ASYMMETRIC PRODUCTION OF THE DM RELIC 
DENSITY 


As said above, in this Section and the following one we 
stick to the minimal IDM framework, as defined in the intro¬ 
duction. Let us make two simple starting assumptions. First, 
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let us assume that the symmetric component of the relic den¬ 
sity left after freeze-out is smaller than the observed value. 
Fast SM gauge scatterings automatically care for that for m^o 
within the ~ 120 — 540 GeV range, whereas for other val¬ 
ues of nifjo large enough A, 3.4 interactions can take care of 
that in. This implies a symmetric annihilation cross sec¬ 
tion larger than the usual thermal freeze-out value ~ 1 pb, 
which means a freeze-out temperature Tfo smaller than the 
usual Tfo ~ nijjo /25 value. Second, let us assume that aB~L 
asymmetry has been generated at a temperature Tb-l above 
m //2 and above the EWPT temperature Tew (which we take 
as the temperature where the vacuum expectation value of the 
SM scalar field becomes sizable, that is Tew ~ 165 GeV from 
Ref 113). We do not care about the way this B — L asymmetry 
could have been generated. It could be due for example to the 
straightforward leptogenesis mechanism. Note that, as well- 
known, if a B — L asymmetry is generated at high temperature, 
a Hi-Hl asymmetry will also be created automatically at high 
temperature from thermal equilibrium SM interactions HI. 

If a B — L (and thus H\ asymmetry) is created at high tem¬ 
perature, an inert doublet H 2 -H 2 asymmetry is to be expected 
too. The scalar potential of Eq. 0 contains the A .5 interaction 
which uniquely does not conserve the number of H 2 minus the 
number of (as well as the number of H\ minus the number 
of HI). This interaction is in thermal equilibrium at T ~ 
if at this temperature the corresponding T^^^ scattering rate, 
given in the Appendix, is larger than the Hubble rate, which 
gives the condition 

l?^5|>10-'’-(m,/,/TeV)'/^ (5) 

If Eq. (|^ is satisfied, the A ,5 interaction equilibrates the H 2 
and H\ (and B — L) asymmetries.* In particular even if, as 
we assume here, no H 2 asymmetry is created at high ener¬ 
gies, such an asymmetry will be created anyway as soon as 
the B — L asymmetry is created. In other words, the inert 
DM model contains an interaction which basically implies 
that “Higgsogenesis” M production of a DM asymmetry is 
at work.^ Note that the scenario could work also the other way 
around, i.e. a primordial DM asymmetry could be at the ori¬ 
gin of baryogenesis via the same X 5 equilibration interaction, 
a possibility we will not consider here (for a scenario of this 
kind see Ref. 113 ). 


' Actually, in the few-TeV asymmetric inert DM scenario considered in 
Ref. 1151 . it is assumed instead that the ^5 interaction could have never been 
in thermal equilibrium, in this case, the DM asymmetry would have been 
created explicitly at high energies, basically independently of the B — L 
asymmetry, if the inert scalar is the DM particle, it turns out that the lower 
bound of Eq. implies that Eq. must anyway hold for TeV masses. For 
instance, Eq. Q with a 100 keV (180 keV) mass splitting implies Eq. 
for m//, > 100 GeV (30 GeV). 

2 in Ref [H, a Xi fermion singlet DM framework is considered with an 
extra X 2 fermion doublet and an X-symmetry. An X 2 asymmetry is created 
from aX-symmetry violating Xl^/f non-renormalizable interaction, which 
is afterwards reprocessed into a X\ symmetry through X 2 decays. Asym¬ 
metric frameworks based on the equilibration of the SM scalar asymmetry 
with a dark sector asymmetry, based on several new dark sector particles, 
or based on various possibilities of a SU (2)i multiplet, can also be found 
in Refs, tni and ca respectively. 


In the following, we will consider in details and chrono¬ 
logically what happens when the temperature of the Uni¬ 
verse cools down from T ^ iuei^ to today T ^ Tew^ crossing 
OT //2 > Tij > Tfo > Tew, with Tx^ the temperature where the 
scattering induced by the X 5 interaction decouples and Tfo the 
freeze-out temperature at which the total annihilation cross 
section decouples. Given that the inert doublet components 
undergo gauge interactions, Tf^ is sizably larger than Tfo, un¬ 
less A-s is of order one, which as we will see is not a viable 
option for the case we are interested in (where the DM asym¬ 
metry is responsible for most of the relic density). Similarly, 
as we will see, Tfo > Tew, i-S- few-TeV DM, is also gener- 
ically necessary in order to have a viable scenario (as in the 
scenario of Ref. ifTSll l. but some violation of this inequality is 
possible. A sketch of the scenario, applied to our framework, 
is shown in Eig.[T] 


A. T > inH2 

At temperature above Tew, 4 inert doublet compo¬ 
nents have a common mass = m 2 . If Eq. Q is satis¬ 
fied, the chemical potential of both scalar doublets are equal, 
BH 2 — BHi ■ Together with the usual SM chemical equilibrium 
relations (from thermal equilibrium SM processes iflTll l. the 
= p .^0 relation (from e.g. r|+r|*** o SM processes), and 
the hypercharge relation 

H (^2; + 4Ah, — 2Adi — Ai- — 2Ae,.) -f Ahi + A//^ = 0 , (6) 

i 

it simply gives 

4 

A/Tj = A//j = ■ (7) 

Erom now on, we define for each species X the asymme¬ 
try Ax = Yx — Yx and the total density 'Lx = Yx +Yx, where 
Yx = nxIs is the particle number density-to-entropy ratio of 
X. Since we are dealing with asymmetries, we also define the 
number of degrees of freedom by summing the number of par¬ 
ticles (or antiparticles but not both), i.e. gx = i for a SU{2)e 
singlet, and gx =2 for a doublet. 

As well-known, for similar B — L and DM asymmetries, 
the DM relic density constraint requires iueim to have a mass 
of few GeV (more exactly, from Eq. 0 and taking into ac¬ 
count the Yb-l to Yb ratio which holds in this case, Eq. Q 
below, one would need mE)M ~ 10 GeV). As this possibility 
is excluded by collider constraints, this implies that a sub¬ 
sequent suppression of the DM asymmetry by a factor of 
~ (10 GeV/moM) must necessarily occur. Two different 
types of suppressions can naturally take place. A first one is 
a Boltzmann suppression from asymmetry violating scatter¬ 
ings, used in several other DM models, see e.g. Refs. lfT^l20l . 
In our scenario, it can arise from the A 5 interaction within the 
period mH 2 > T > Tx^. The other possible suppression can 
arise later when T < Tew from the combined effect of DM 
oscillations and symmetric annihilations. 
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B. 


mHz ^ T > 


Once the temperature drops below m/Zj, if the X 5 interaction 
goes on to be in thermal equilibrium, the H 2 asymmetry gets 
Boltzmann suppressed. This can be directly seen from the 
Boltzmann equation of A// 2 , valid for T > Tew, 


dAH 2 

4 

/A//2 

A//i 

dz 

sHz 


ye? 


( 8 ) 


where z = mH 2 /T, H(z) is the Hubble rate and is the 

reaction density of the Xs scatterings, given in Appendix. It 
includes both pair annihilation/creation H 2 H 2 O H\H\ and “t- 
channel” H 2 H 1 O H\H 2 processes. The interaction leaves 
intact the sum of the asymmetries of H\ and H 2 but not each 
asymmetry individually. Once T drops below the hrst 
term in the r.h.s. of Eq. (| 8 ]l is enhanced with respect to the 
second term by the fact that is Boltzmann suppressed, un¬ 
like This Boltzmann suppression of the asymmetry lasts 
until the Xs induced scatterings decouple, at T = T-^^, when 
~ Quantitatively, this can be accounted by the usual 
k-factor which gives the asymmetry as a function of the tem¬ 
perature 


Ah2 = gH2 A'//2 , (9) 



Down to , the chemical potential relation /J //2 = still 
holds, and we get 


^H 2 {z) 



( 11 ) 


This is nothing but the solution which makes the r.h.s. of 
Eq. ([^ to vanish. Together with the other chemical potential 
relations above, the asymmetry reads at T = (similarly to 
the fermion doublet case of Ref. 061) 


^H 2 {zx^) 


16k 

158+13k (z^,) 


( 12 ) 


Eor practical reasons, it is convenient to dehne 

a| = |Ah 2 (z^,)|. (13) 

Let us note that since ^ 1, the k-factor can be approxi¬ 
mated by 


(14) 


Clearly, the Xs coupling must not be too large in order to avoid 
a too strong exponential suppression of the H 2 asymmetry. 

As emphasized in Ref. DU, for fermion quartic interac¬ 
tions, the last A .5 induced channels to decouple are the t- 
channel ones,^ simply because they are less Boltzmann sup¬ 
pressed than the other ones. The value of z;^.^ is given by 


the condition that the E^^^ rate is equal to the Hubble rate H. 
One has for example that if m//; = 10 TeV, z^j = 2(15) for 
A ,5 = 3 • 10^®(5 • 10^®). Eor somewhat smaller values of A, 5 , 
even if the A .5 reactions doesn’t enter in thermal equilibrium 
(as dehned by Eq. 0 ), a numerical integration of the Boltz¬ 
mann equations shows that still a number of scattering pro¬ 
cesses occur nevertheless, what can lead to a sizable asymme¬ 
try. 


C. > T > Tfo 

During this period, the Ah 2 asymmetry stays constant un¬ 
like the total abundance E// 2 , whose Boltzmann equation for 
this period reads 


d'LH2 ^ {^effv)s r 2 _^5Z_^eq2 
dz ZH V "2 H 2 H 2 


(15) 


where {Oeffv) is the effective thermal cross section of the 
H 2 H 2 ^ SMSM annihilations, given in the Appendix. With 


a constant A^^, as it is the case during this period, the solution 
of Eq. (15 1 at freeze-out is to a good approximation given by 


^^ 2 (^/ 0 ) — 



2 

-H 2 



(16) 


which is nothing but the expression which makes the r.h.s. of 
Eq. ( |T5| ) to vanish. Here, by z/o we mean the usual freeze-out 
value given by the equation 


0.003S-mpi2gH2mH2{<yeffv) 

2/0 - In-— -- 

y/gTzf^ 


(17) 


If the annihilations are fast enough to leave at Tfo a sym¬ 
metric component smaller than the asymmetric one (which 
is typically satisfied for (Oe//v) > 1 pb), the following rela¬ 
tion holds, Eh 2 (z/o) — ^2 ^ ^H 2 (^fo)- Given the sign of the 
baryon asymmetry, this means at Tfo, E //2 ~ —^H 2 ~ -Yh 2 > 

Yh 2 XI- 


D. Tfo > T > Tew 

Nothing is expected to happen during this period. The H 2 
total density left at Tfo is left intact until Tew, temperature 
at which the total density and asymmetry are given by (for a 
dominant asymmetric component) 

Y.H2 {ZEW ) — S//2 (2/0 ) — A^2 (ZEW ) I = A^2 • 

(18) 


E. T < Tew 


2 We thank G. Servant and S. Tulin for discussions on the importance of Next, once the temperature drops below Tew , tWO new ef- 

these channels. fects enter into play: generation of mass splittings between 
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the H^\ and ri+ components and possibly fast inert particle- 
antiparticle oscillations r]*' r]*'*. The effect of the mass split¬ 

tings generated by the SM scalar vev, Eq. Q, is of moderate 
importance. Assuming, as said above, that the component 
is the lightest one (i.e. < 0 ), they imply that the other com¬ 

ponents will ultimately decay to H^. But these decays con¬ 
serve the number of inert scalar particles. They just convert 
the H 2 asymmetry created before EWSB (with mass m/z^) into 
a DM relic density of selfconjugated DM particles (with 
mass = nioM, different from unless Xjjo vanishes). 
More important is the potential effect of the much faster in¬ 
ert particle-antiparticle oscillations r|‘' o r|®* caused by the 
A ,5 interactions. The rate of an oscillating particle is simply 
given by the value of the associated mass splitting osiiniEa, 
i.e. Toic = 5m = m^o — m^o. Eor T < Tew^ this rate is very fast 
compared to the Hubble rate 


H{T) 


~ 2 - 10 '^ 



lOOGeV 

f 


\mHoJ 


(19) 


The effect of these oscillations depends obviously on whether 
they do occur, which Eq. ( [T9| doesn’t necessarily imply, and 
on whether symmetric annihilations do occur after EWSB. 
Actually, even if the freeze-out occurs before EWSB, this does 
not imply that symmetric annihilations could not restart again 
after EWSB, due to oscillations Emm. This could easily 
be the case because, even if one starts with a pure asymme¬ 
try, the oscillations will quickly give a number density of each 
population much larger than their thermal equilibrium values, 
roughly n^o ^ ^ | An^o | /2 ^ , so that | An^o | (ov) > H 

can hold even if (ov) < H. If these annihilations occur, 
they will anyway reduce the DM abundance, as no inverse 
processes will occur in this case. Let us consider both possi¬ 
ble cases separately. 


1. T < Tew No symmetric annihilations after EWSB 


If no symmetric annihilations arise after EWSB, oscilla¬ 
tions have simply no effect. They quickly reconvert a pure r|*’ 
population, or a pure population, into an oscillating mixed 
r| 0 -r| 0 * population, but they do not change the number of inert 
states ED- In this case, the number of particles left today 
will be simply equal to the number of inert scalar particles 


stored in the H 2 asymmetry before EWSB, i.e. ~ A^^, 

that means the density is equal to the asymmetry left af¬ 
ter A -5 interaction’s decoupling. Erom Eqs. (12 1 and (18 1 , this 
gives 


ytoday 
^ DM ' 


16k {zx^) 

158 + 13k//2 (^^ 5 ) 


Ab-l , 


( 20 ) 


with k{zx^) given by Eq. (14 1 . Only the relation between the 
value of Yg today and Ab_e changes after EWSB, as a result 
of the fact that below Tew the conservation of electric charge 
holds rather than conservation of Y and Tj. We get 


ytoday 


12 


Ab-l ■ 


( 21 ) 


As a result, the DM density reads 

ytoday _ 148'k(z^j) ytoday 

“ 474-f39-k(z^J ' 


( 22 ) 


and the actual DM to baryon density ratio is given by 

^DM _ / niuo \ _ 148-k(z^J / \ 

~ ’ VI GeV/ “ AlA + 39-k{zXi) ' VI GeV/ ’ 

(23) 


This is the final result if no symmetric annihilations occur af¬ 
ter EWSB. However, it must be stressed that it is not manda¬ 
tory to avoid symmetric annihilations after EWSB. On the 
contrary, if the Xs interaction above does not provide enough 
suppression, these scattering processes could easily provide it, 
without the need of any special tuning. This is what we will 
now quantify. 


2. T < Tew- Symmetric annihilations after EWSB 

Possible effects of dark matter oscillations have been stud¬ 
ied in Refs. ED ED. As said above, since oscillations repro¬ 
cess the asymmetry into both particle and antiparticle densi¬ 
ties, their main effect is to allow the symmetric annihilations 
to start again. Even if, as Eq. ( [T9| ) shows, the oscillation rate 
is much faster than the Hubble rate, this doesn’t necessarily 
mean that oscillations (and thus eventually annihilations) do 
occur. As shown in these references, if dark matter under¬ 
goes fast annihilations or elastic scatterings, these processes 
can break the coherence of the r|®-ri*** states, preventing them 
from oscillating. The interplay of the oscillations with the 
other processes is actually more complicated in our doublet 
scenario than in the singlet setups considered in Refs. Eni22l . 
On top of annihilations and r]*' or ri*** elastic scatter¬ 

ings, there are charged r|^ states, which at Tew are respon¬ 
sible for half of the asymmetry and do not oscillate. East in¬ 
elastic scatterings can change neutral states into these charged 
states and vice et versa. Moreover, as said above, all states 
ultimately become real states, which obviously do not os¬ 
cillate. 

Let us first consider what happens to the neutral states, as 
if there were no charged states. In this case, one has two 
important processes. On the one hand, there are 
SMSM annihilations processes which are dominated by their 
A. 3_4 interaction contribution. On the other hand there are 
r|0(*)5M — 7 > elastic scatterings which are dominated 

by their t-channel Z exchange contribution. The later domi¬ 
nates over the A ,3 4 elastic scattering contribution. Eor T > m/,, 
this stems from the fact that it involves a t-channel media¬ 
tor whose mass is much smaller than the inert scalar mass, 
mz ^ m/Zj. It scales as ~ G\T^ as compared to the 

quartic coupling elastic contribution which scales as r^cat"‘^ = 

^ ^ *- 6 ® quartic contribution 

is also subleading because it is Boltzmann suppressed, unlike 
the gauge one. As a result, the gauge elastic contribution is the 
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last to decouple. The Z exchange process is relevant for pre¬ 
venting states from oscillating because the gauge inter¬ 
action is odd under exchange ll^ . Thus the relevant 

question is, down to which temperature will these processes 
effectively prevent the oscillations to start ? At first sight, we 
could think that oscillations will start only once the scatter¬ 
ing rate Tscat goes below the oscillation rate To^c = 5'^* (if 
this time both rates are still larger than the Hubble rate). This 
turns out to occur at a rather low temperature, Tqsc few GeV 
scale. In this case one would be back to the “no symmet¬ 
ric annihilation” case above, because oscillations have practi¬ 
cally no more effect at this temperature, where the annihila¬ 
tion rate is already largely suppressed. However, an integra¬ 
tion of the Boltzmann equations shows that oscillations rather 
start when = Tscat, see Ref. ETIl . In our scenario, 

as we also have checked from a numerical integration of the 
relevant Boltzmann equations, this turns out to happen at a 
temperature above Tew- Thus we conclude that oscillations 
start as soon as EWSB occurs. As a result, annihilations can 
restart from this temperature and to determine how much of 
them will annihilate, one can just take the Boltzmann equa¬ 
tions with the oscillation and annihilation terms. 


dz 

c/A^o 

dz 

£/S^o 

dz 


S ry 2 
2zH I 


0 



yeq 2 

S® J ’ 


.5m „ 

= 2 ;— ^^0 , 
-.5m 

= 2;— A^o- 


(oov)i „ 


(24) 

(25) 

(26) 


where for any T < Tew we define z = m^o/T, with (cov) 
the thermally averaged —>^ SM annihilation cross sec¬ 

tion, and a quantity that accounts for the coherence be¬ 
tween the T]*' and 'q*'* components (see ED for further de¬ 
tails). The resolution of these equations leads to a monoton- 
ically decreasing E^o (z) function and to oscillating functions 
A^o(z) cos[/(z)] and S^o(z) sin[/'(z)] whose amplitudes 
also decrease monotonic ally. For fast oscillations, and ne¬ 
glecting the term in Eq. (24 1 , the set of Boltzmann equa¬ 
tions can be simplified and solved analytically, at an approxi¬ 
mate level, as explained in the Appendix. The solution it gives 
for E^o is ^ 


2 

Eno (z > zew) = -——TTT 

' I I 1 (gov)^(z) 




r97'l 


with Zew = /Tew and where we fixed the initial abun- 

dance and asymmetry to be equal to E^o(z£w) = ^//t/ 2- The 
asymmetry A.p 0 and S.^o are, in turn, fast oscillatory functions 
which are equal to zero on average. 


^ This result is approximately the same than the one obtained in 12 T1 for 
much smaller Sni values — see Eqs (25) and (33) therein - but in which 
Xcsc,ann (whlch depends on Sm) is now simply replaced by ■ 


The result of Eq. ( p7| can also be qualitatively understood 
in the following way. Once T < Tew, the fast oscillations re¬ 
process quasi instantaneously the asymmetry in oscillatory 
abundances for q*' and q***. On average, just after EWSB, we 
have therefore n^o ~ n^o, ~ lAn^olr^^/2. Since Tfg > Tew, 
when two conjugate particles annihilate to SM particles, the 
reduction of inert doublet state it implies will not be compen¬ 
sated by any inverse processes. As a result, the Boltzmann 
equation for E^o one gets along this way is simply given by 


c/E^o 

dz 


y 2 

2zH ^'1“ 


(28) 


whose resolution leads to nothing else than Eq. (27 1 .^ 

The next step is to include the contribution of the charged 
states. Since these states do not oscillate, one could naively 
expect that the charged asymmetry is essentially left intact 
until the charged states decay to states. This doesn’t 
work this way. To see that precisely, one should in prin¬ 
ciple solve the corresponding set of six coupled Boltzmann 
equations, for Et^^, A^^, Hrm, Et^^, At^^, where Hi^g and 
Erio are the real and imaginary parts of the quantity that ac¬ 
counts for the coherence effects. Nevertheless in practice 
we don’t need to go that far. It turns out that, if just before 
EWSB there are essentially only q^ and q*’* states as con¬ 
sidered here, as soon as oscillations start they put the neutral 
state asymmetry to zero (on average), and processes which 
can transfer a charged asymmetry into a neutral one will very 
quickly put the charged asymmetry to zero too. This will be 
done in particular by q*’(*) SM q^ SM inelastic scatterings 
and q* o r\^^*^SMSM decays. The decrease of the charged 
component asymmetry due to these processes is exponential 
(A^+ oc e as can be seen from the corresponding term in 
the Boltzmann equation, sHzdS.^+ /dz —A^^+y^+^^o -I-...). 
This “re-equilibration” of the asymmetries by these processes, 
which follows their “desilagnement” by the oscillations when 
these latter start, occurs much faster than the process of sup¬ 


pression of Eria in Eq. (50 1 . Asa result, in the same way as for 


the neutral states, one can adopt the simple assumption that as 
soon as oscillations start, the particle and antiparticle densities 
for charged states are equilibrated, = V=V. At 

this point, the annihilation processes such as q+q^ —> SMSM, 
q+qO* g]\^g]\^ and q^q*^ —> SMSM can start again, in the 
same way as the q^’q*'* —SMSM ones. The whole effect can 
be approximatively accounted by the simple Boltzmann equa- 


^ The reason why the two results coincide is in fact more subtle. Since the 
oscillatory behavior is given by 

the Boltzmann equation in this naive approach should read 


dz 



1^0 TpO* 


(gpr)^ 2 

2zH 


sin^g(z) . 


Averaging this expression, we find Eq. j28| up to an extra factor 1/2. An 
extra factor 2 must nevertheless be added to take into account the contribu¬ 
tion of the coherence part, giving back Eq. 0- 
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d^H2 _ {^effv)s ^2 

~dr~ 


(29) 


Similarly to what has been obtained in Eq. (27 1 , the resolution 
of Eq. ( |29l l, integrated from Tew until now and using the initial 
condition in ([TSl), leads to 


S//2(2 > ^Ew) 


X 

^H 2 


1 , {<^effv)s{z) 

^ + H{z) 




(30) 


This equation holds for the case where the total number den¬ 
sity just before EWSB is given by the asymmetry. If there 
is also a non-negligible part which is left from the symmetric 
freeze-out, one must simply replace the asymmetry at zewi 
by the total number density at the same temperature, 
^H 2 since this is the number which determines the num¬ 

ber of symmetric annihilations which will occur after EWSB,® 


'^H2 {z > Z£w) 


1 -f 


^H 2 {zew) 

(o^//v)»(z) f ^ _A 

H{z) J 


^H 2 Uew) 


(31) 



ZAs Zfo ZEW 


FIG. 1: Sketch of the scenario considered in section II and III. We 
represent, as a function of z = lUfjo/T, the H 2 asymmetry | in red 
and the total DM density in black. We neglect in this sketch the 
mass splittings between the different components of H 2 - First step: 
The initial asymmetry and A//^ are fixed by the B — L asym¬ 
metry. Second step: this asymmetry gets suppressed until the 7,5 
interactions decouple at zx^ — the suppression is characterized by 
the fc-factor. Third step: at zx^ the total annihilation cross section is 
still in thermal equilibrium and follows the equilibrium density 
until it reaches |A/fj|. Fourth step: at EWSB, oscillations start 
and reequilibrate the particle-antiparticle populations. At this point, 
annihilations can start again and if they do they deplete the density - 
the suppression is characterized by the K-factor. 


with E// 2 (z£w) — '^H 2 {zfo) as given in Eq. ([16 1 . 

Note also that the rfSM o r[^ SM (and conjugated) pro¬ 
cesses above not only equilibrate the neutral and charged 
asymmetries, but also can break the coherence of the 
by transforming a coherent neutral state into a charged state 
which does not oscillate. However, in the same way as for the 
Z exchange channel above, its rate goes under {dm)^/H be¬ 
fore EWSB occurs, so that they do not prevent oscillations to 
start at r = Tew ■ 


Einally, because of the mass splittings between and the 
other components of the inert doublet, this total density is pro¬ 
gressively transferred into a H'^ density through the decays of 
the heavier components. Note that, as we will see, in the nu¬ 
merical section below, Eq. ( (3T] i reaches its asymptotic value 
to a good approximation before T drops below the value of 
the mass splitting — nifjo. As a result, this splitting can be 
neglected as it was done to get Eq. 


F. Final inert scalar relic density 

We summarize in Eig. [T] the evolution of the asymmetry 
IA// 2 I and total density We remind the main steps: 


® If there is no asymmetry and if the freeze-out has occurred prior to EWSB, 
one recognizes in Eq. the usual asymptotic freeze-out behavior, i.e. the 
freeze-out is not instantaneous, but reaches asymptotically its final value as 
given in this equation. In practice, as well known, the effect is negligible 
in this case, i.e. the denominator is equal to unity to a good approxima¬ 
tion. Here, instead, the denominator at Zfo can be much larger due to the 
asymmetry. 


l. T > niH^. The H 2 asymmetry, proportional to the B — L 
asymmetry, is generated through the A-s interactions: 

^H2{z<1) = -^Ab-l. 0 

2- >Tx^. The asymmetry undergoes a Boltzmann 

suppression until the A 5 interaction decouples, 

K = (^^ 5 ) I = + ® 


3. Tx^> T > Tew- The symmetric component of E //2 fol¬ 
lows an exponential suppression, until reaches at 
z = Zfo the value 


^H 2 {Zfo) - 



peq 2 
-H 2 



ini' 


Erom Zfo to Zew, the annihilations are momentarily 
frozen, so that T,h 2 {zew) — ^H 2 {zfo)- During this pe¬ 
riod, if the contribution from usual freeze-out is neg¬ 
ligible, there are only H 2 particles in the plasma and 

^H2{zew) — A^j. 

4. T < Tew- The fast r|‘’ o r]®* oscillations start. They 
quickly reprocess the H 2 density in equal abundances 
(on average) for r]*', r|®*, r|+ and r|^. The annihilations 
can therefore start again and, if they do, they deplete the 
set of densities, whose sum reads asymptotically 


E//2(z » zew) = ' 


'^H 2 {zew) 


1 , -liz) c V ^ 

1+ H{z) I^^H 2 [Zew) 


(32) 
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This total density is progressively transferred into a 
density through the decay of the heavier components. 


The final DM abundance is therefore given by 


^today V i'-^ 


where we define 

{^effv)s{z) Z 


K = 


H (z) zew 


~ 1.3-10 


13 


\^effv) 

1 pb 


(33) 


(34) 


Since ultimately no asymmetry survives, the relation between 
the baryon and the B~L asymmetry is still given by Eq. ( [2l] i, 
and the final DM to baryon density ratio is given by 

^DM _ ^H2{zew) 1 / 'W/fO \ 

Hb “ l+K-'LH 2 izEw) ' ■ U GeVJ ^ 


or equivalently, if the asymmetric component dominates, us¬ 
ing Eqs. ( [T2 ) i and ( [T8] l, 


Q.DM 

Q.B 


_ 148 _ 

474 + (39 + 148kTb"''"^) k {zx,) 



A number of comments can be done regarding these results; 


Eqs. (33 I and (35 1 show that beside the Xs interaction 
induced “A:-factor” suppression in A//^, see Eq. ( [T^ , os¬ 
cillations drive a 1 -\-k-'Lh 2 {zew) factor suppression. 
This “K-factor” suppression can be sizable as soon as 
^■'^H2{zew) ^ 1 - 


• As Eq. shows, this suppression is neither instan¬ 
taneous nor exponential. It goes as the inverse of 
zjzEW ~ 1 until it reaches an asymptotic value. In this 
sense, imposing that the cross section satisfies the uni- 
tarity bound, it is naturally limited but it still can be re¬ 
sponsible for the ~ {IQGeV/mEiM) suppression needed, 
see below. 


• The appearance of the K • (zew ) factor is not surpris¬ 
ing. The condition K • T,h 2 {zew) < 1 is nothing but the 
condition (n//^ + ^H2){'^effv) <H aXT = Tew- 

• Interestingly, for large values of K • E/Zj {zew), the ^dm'^ 
relic density obtained doesn’t depend anymore on the 
asymmetry left at Tew, even if this asymmetry is the 
source of the final DM abundance. In this case, we sim¬ 
ply get 


ytoday _ ^ 

^DM ' ~ ’ 


and the ratio reads 


Q-dm 

TIb 


0.15 -zbw ■ 



(37) 


This means, as we could have anticipated, that for large 
cross section the asymmetry left is independent of the 
initial asymmetry, provided this initial asymmetry is 
large enough. In other terms, if the K-factor suppression 
is small, both baryon and DM asymmetries are directly 
connected. If instead it is large, they are not related any¬ 
more in a so direct way, since in this case the final relic 
density depends only on the annihilation cross section.^ 
Note interestingly that Eq. ( [37] i is nothing but the result 
of the standard freeze-out scenario, but with the impor¬ 
tant difference that in the standard case, zew in Eq. ( p7| ) 
must be replaced by z/o. 


III. FAILURE OF THE ASYMMETRIC IDM SCENARIO 


The final result of Eq. ( p5| ) depends on three parameters; 
m^o, E// 2 (z£w) and the total cross section (Ce//v) via K in 
Eq. ( [34] l. This means that for given values of the input param¬ 
eters nifjo and {Oeffv), there is only one value of '^h 2 {zew) 
which gives the observed value of TIom/^b, as given by the 
PLANCK best fits, = 0.120 and TIbH^ = 0.022 H]. 

Since T.H 2 (zew) depends only on these two input parameters 
and on this means also that there is only one value of 
A^^ which gives the correct relic density for fixed values of 

the two input parameters. We show in Eig. 2 this value of aJ^^ 
as a function of m^o for different values oHhe cross section. 


1 

By comparing this value of A^^ to the value this asymmetry 
would have if there were no “A:-factor” suppression - given by 
the A//j upper horizontal line - one can read off what is the 
value of this A .5 induced “A:-factor” suppression, Eq. (12i as 
compared to Eq. 0 . 

As said above, to dominate the final relic density, the 
asymmetry cannot be suppressed by more than a factor 
'«dm/ 10 TeV. Eigure|^also shows the corresponding values 
of the K • E //2 {zew) factor which lead to the other suppression, 
i.e. the 1/(1 + K• E/z^Czew)) factor in Eq. (35i. It also shows 
for which values of the various parameters the asymmetry pro¬ 
duced before the EW transition is responsible for 50% of the 
final DM relic density (black line). Above (below) this line the 
relic density is dominantly of asymmetric (symmetric) origin. 
Similarly, the dotted upper (lower) black line gives the val¬ 
ues of the parameters above (below) which the asymmetry is 
responsible for more (less) than 90% (10%) of the final relic 
density. Eor masses which give a freeze-out below Tew, the 
k-E// 2 (z£w) factor becomes exponentially large because in 
this case T,h 2 {zew) is still exponentially larger than its value 
at freeze-out. Thus, the proportion of T-fj^izEw) which is due 


to aJ/j is therefore exponentially suppressed. This explains 
why the black lines quickly go up for moM below 4 — 5 TeV. 
Note nevertheless that this suppression, even if exponential, 
is far from instantaneous. As a result we find that, still, the 


^ But still, even in this case, they remain similar as the K factor is bounded 
from above by unitarity considerations on the total cross section. 
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function of the input parameters m^o for various values of {Ogffv) 
(dashed blue lines). The corresponding values of k-£// 2 (z£w) ^6 
given by the dashed orange lines. The upper horizontal line gives 
the value of which is obtained from equilibration with the Hi 
and B — L asymmetries. The r.h.s. solid (dashed) green line gives the 
maximum value of the input parameters imposing that A .3 4 couplings 
are smaller than An (\/ 47 t). Below lUfjo ~ 4.7 TeV, the freeze-out 
occurs after EWSB, which relatively quickly causes a huge suppres¬ 
sion. The black lines from top to bottom give the value of the param¬ 
eters for which 90%, 50%, 10% of the relic density is of asymmetric 
origin, respectively. 


asymmetry can dominate the relic density for a mass equal to 
3.7 TeV which is substantially lower than the 4.7 TeV value 
which gives Tfo = Tew-^ A comment which must be made at 
this point concerns the fact that we have considered the elec- 
troweak phase transition as if it was an instantaneous process, 
i.e. as a step function at the temperature Tew ^ 165 GeV - 
from Ref. ||T3 (see also Ref. M) - which as said above is 
the temperature where the vacuum expectation value of the 
SM scalar field becomes sizable (i.e. where the oscillations are 
about to start to reprocess the asymmetry). As the electroweak 
transition is a crossover, it is clearly an approximation which 
could be refined. A change of Tew by a given factor would 
shift all values in Fig.|^by about the same factor. 

As expected from the discussion above. Fig. also shows 
that, for large value of K • { zew )^ the observed relic density 
doesn’t depend anymore on the value of A^^, provided this 
later quantity is above a certain value. 

Note that the r.h.s. green curve of Fig.|^is obtained by im¬ 
posing that all quartic couplings are perturbative, A ,3 4 < 47t. 
This line shows that a dominant asymmetric component re¬ 
quires that nioM ^ 25 TeV (whereas the same condition gives 
tnoM 30 TeV for the standard freeze-out scenario and for 
a small value of the coupling). Such a bound also im- 


* To get this 3.7 TeV value we simply applied Eq. j30| neglecting the fact 
that in this case the inverse scattering term must be taken into account 
in the Boltzmann equations (as in Eq. 124) ). The incorporation of this term 
would slightly lower further this minimum value of nioM- 



10 

1 

10 '^ 

10'2 

10 ’^ 


> 

o 




1 

FIG. 3: Values of which leads to the values needed in 
Fig. HI as a function of the input parameters m^o. The red shaded 
area is excluded by the direct detection constraint of Fq. taking 
m^o — m []0 = 100 keV. Taking instead m^o — m^o = 180 keV gives 
the red dashed line. The dashed black line gives the value of 7,5 be¬ 
low which this interaction never gets in thermal equilibrium, as given 
by Fq. 0. The lowest allowed 7,5 value is obtained for a situation 
where there is neither a k suppression, nor a K suppression. In this 
case one generates directly the observed relic density from having 
only partial thermalization of the asymmetries. 


plies an upper bound on the (Oe/fv) cross section of about 
2.5 pb, that is to say a value about 4 times larger than the 
~ 0.7 pb value one needs at these energies along the standard 
freeze-out scenario. Imposing instead that 7,3 4 < one 
gets moM < 8 TeV and (Oe//v) < 1.1 pb (dashed green line). 

The minimum value of the Xj -I- 7^ coupling combina¬ 
tion (which enters in (Oe//v)) that this scenario requires is 
~ 2, corresponding to iiieim ^ 4 TeV and a cross section of 

~ 0.5 pb. This is smaller than the usual ~ 0.7 pb because 

1 

the associated asymmetry ^h\ ~ |A//j I also participates to the 
depletion of the total density. No need to say that with such 
large values of these quartic couplings, Landau poles are to 
be typically expected far below the Planck scale. Although 
the energy scale at which we get a Landau pole depends on 
the value of other couplings such as 72 , if there is no can¬ 
cellations between the contributions of various couplings in 
the beta functions, a value of 73^4 ~ 1.5 gives a Landau pole 
at ~ 10^-10® GeV. This means that new physics is to be ex¬ 
pected in this case below this value. The scale of B — L asym¬ 
metry production has not to be necessarily below this scale. 
All what matters for the value of Q.dm /is the value of the 
B — L asymmetry at T ~ 111^2 ■ 

In Fig. 0 as a function of the same two input parameters 
m^o and (Oe//v), we show the value of X 5 which leads to the 

A^^ value needed in Fig. The corresponding value of the 
mass splitting is also given on Fig. 0 This figure 
shows that the scenario leads to the observed renc density for 
75 S [5 • 10^^, 8 • 10^®], which corresponds to a mass splitting 
equal to approximately ifiAg — mhq S [0.1,15] keV. Larger val¬ 
ues of 75 quickly lead to a X^ decoupling temperature much 
smaller than thus to largely Boltzmann suppressed re- 
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FIG. 4: Same as Fig.|^ but allowing the inert scalar to decay into a 
lighter real scalar singlet S with mass = m^o/10. 


FIG. 5: Same as Fig.|^ but allowing the inert scalar to decay into a 
lighter real scalar singlet S with mass = m^o/10. 


maining asymmetries. Smaller values of A .5 rather quickly 
lead to no thermalization of the H 2 and Hi asymmetries, i.e. to 
no creation of a H 2 asymmetry. In most of the relic density al¬ 
lowed parameter space, both the “k” and “k” suppressions are 
active, although it is possible to have only one of the effect to 
account for all the necessary suppression. As said above, an 
important constraint that one must satisfy is the direct detec¬ 
tion constraint of Eq. 0 . The value of the mass splitting just 
quoted are below the ^100 keV direct detection lower bound 
of Eq. 0. 

Thus, unless direct detection would allow a mass splitting 
as low as the value m^o —irijjo ~ 15keV, which seems very 
unlikely, this very minimal asymmetric scenario is in fact ex¬ 
cluded ! This can also be clearly seen from Eig. 0 where the 
region allowed by direct detection taking in Eq. 0 a mass 
splitting dirimin = 100 keV has no overlap with the region 
which gives the observed relic density. Or, in other words, 
imposing that the mass splitting is above lOOkeV, the in¬ 
teraction turns out to decouple only at zx^ ^ 50 leading to a 
tiny Q.DM relic density. 


IV. REPROCESSING THE INERT DOUBLET 
ASYMMETRY INTO A LIGHTER PARTICLE DM RELIC 
DENSITY 

Since the very minimal IDM scenario above cannot account 
for both the relic density and direct detection constraints at the 
same time, one question one must ask is whether this simple 
scenario of an IDM asymmetry creation could not be never¬ 
theless at the origin of the DM relic density today in a simple 
way. This could in fact happen if the DM is made of a lighter 
specie, whose relic density would be due to the reprocessing 
of the inert doublet asymmetry into this specie. Such a repro¬ 
cessing could for instance take place through decay. Eor the 
scalar scenario we consider, this could be the case if there ex¬ 
ists a lighter Z 2 odd particle, “5”, to which the inert doublet 
states can decay. In this case, if the asymmetry is fully re¬ 
processed into this lighter particle S, so that each inert scalar 


component gives one S particle, the results of Eigs.|^and[^ 
are still fully valid provided the mass of the S particle, ms, 
is close to m/Zj. If instead it is sizably smaller, this requires 
to create more inert particles by a factor lUfj^/ms- As an ex¬ 
ample, Figs. 1^ and 1^ show the value of parameters we need 
to get the observed relic density for a ratio equal to 

10. Note that such large H 2 asymmetry cases, beside allowing 
smaller DM masses, also give relaxed lower bounds on the 
A, 3_4 couplings (in order to suppress sufficiently the symmetric 
part). Sizably smaller values of these couplings are possible, 
relaxing accordingly the Landau pole constraints. In order to 
reprocess the inert doublet asymmetry into such a S specie, 
various possibilities could be considered. 

A simple possibility is to consider 5 as a Z 2 -odd scalar sin¬ 
glet, into which the scalar doublet states decay sufficiently 
slowly to happen after the freeze-out of this singlet DM par¬ 
ticle. Such a decay can be accounted for by a A 9 ^sHiH 2 S 
renormalizable interaction. If so, the main constraint to sat¬ 
isfy along such a scenario is, in order that the S relic density is 
mainly produced from the IDM asymmetry, that the S particles 
has a S ^ SMSM annihilation channel with a large enough 
cross section to leave a relic density smaller than the observed 
one at S freeze-out. These annihilations can be accounted for 
by a A 9 XsH\HiS^S interaction. 

As an example, if we take a real scalar singlet, and fix 
the parameters to be ms ~ 2 TeV (400 GeV) and m/Zj = 
10 TeV, both conditions are fulfilled for Xs > 0.6(0.1) and 
/Z 5 < 4 • 10^^ GeV (7 • 10^® GeV). Also, the Xs interaction 
induces elastic scattering on nucleon through SM scalar ex¬ 
change, Oyv = where mf] is the nucleon 

mass and the nucleon form factor is approximately given 
by fs! « 0.3. The LUX experiment constraint ll28ll . which 
for moM ^ 100 GeV is Ozz < 1.2 • 10 ^"(m 5 /l GeV) pb, is 
satisfied if < 1.6 • 10 ^^(ot 5/1 GeV)^/^. Combining both 
these lower and upper bounds on Xs leads to the lower bound 
ms > 300 GeV ll29l . In Fig. [^we show the evolution of the 
asymmetries we get as a function of the temperature for an 
example of parameter set which leads to the observed relic 
density. 
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FIG. 6 : Evolution of the various abundances as a function of z = 
m/fj/r in the case where the inert doublet decays into a real state 
S after 5 freezeout. Such an evolution has been obtained fixing 
the H 2 -rel?Lted parameters to m //2 = 10 TeV, (Ogffv) = 0.5 pb and 
5 m //2 = 3 • 10^® GeV (corresponding to X 5 ^ 10^®), the S-related 
parameters to = 1 TeV and (CJ 5 V) = 4 pb, and the connector pa¬ 
rameter controlling the decay rate to = 5 • 10 ^® GeV. The dashed 
purple curve shows what would be the evolution of the 5 density if 
there were no H 2 ^ H\S decay. 


FIG. 7: Evolution of the various abundances as a function of z = 
m/fj/T in the case where the inert doublet decays into a complex state 
S before S freezeout, for m [{2 = 10 TeV, (Ogffv) = 1 pb, 8171^2 = 3 • 
10^® GeV (corresponding to ^,5 = 10^®), = 1 TeV, (o^v) = 4 pb, 

8 ms = 10^^ eV, and /jg = 5 ■ 10^^ GeV. The dashed purple curve 
shows what would be the evolution of the S density if there were no 
H 2 -^ Hi S decay. 


V. SUMMARY 


Similarly to the fermion scenario considered in Ref. 
another possibility is to consider instead that the decays occur 
when the freeze-out of the singlet particle S has still not taken 
place. In this case, the inert doublet asymmetry could also be 
at the origin of the DM relic density, if the singlet is a complex 
field and if the inert doublet asymmetry is reprocessed into a S 
asymmetry. Since inert doublet oscillations start at Tew^ this 
requires the reprocessing to be done prior to EWSB. Imposing 
in addition for simplicity that the decay occurs after the A .5 in¬ 
teraction has decoupled at zx^, for example for mH 2 = 10 TeV 
and ms ~ 2 TeV, one needs 10^^ GeV ^ < 10^^ GeV. For 

this scenario to work, one has to make sure that the S asym¬ 
metry created in this way is not washed-out by possible S-S^ 
oscillations. This requires that terms as Hi {S^ + h.c.) or 
+h.c. are sufficiently suppressed for the oscillations not 
to occur before S freeze-out. This means the S mass splitting 
dms = {mf + Q{Tew — T)'h!sV^)/2ms must be smaller than 


^ 10 ^ 



(38) 


with z^^g > 20 the value of ms/T at which the S freeze-out 
occurs, and (a^v) the S annihilation cross section. Note that 
at temperature lower than when the S oscillations starts 
at they can allow the S annihilation to restart in the same 
way as for the inert doublet above. Similarly to Eq. ( [33| , this 
causes a suppression of the S asymmetry by a factor equal to 
(1 with K 5 = {osv)sz/H{z)zlsc- In Fig. we 

show an example of evolution of the H 2 and S asymmMries 
along such a scenario. 


In summary, if there exists an inert scalar doublet H 2 , unless 
the A-s interaction is tiny, the inert doublet components will au¬ 
tomatically develop an asymmetry from thermalization with 
the ordinary SM scalar doublet and lepton asymmetries. We 
have studied in details what is the fate of such an asymmetry at 
temperature below the value of the inert doublet scalar mass, 
m/Zj. Beside being responsible for the asymmetry creation, the 
A 5 interaction also controls the neutral component mass split¬ 
ting (hence the Z-exchange direct detection rate) and induces a 
“k-factor” suppression of the inert doublet asymmetry at tem¬ 
perature below m/Zj. On top of this suppression one can also 
have an extra “K-factor” suppression, from the combined ef¬ 
fect of DM oscillation (also induced by the A 5 interaction) and 
DM symmetric annihilation. This leads to a scenario which 
chronologically occurs as represented in Fig. We showed 
that in the few-TeV range, there is an all region of parameter 
space where the DM asymmetry survives enough to lead to 
the observed DM relic density, Fig.|^ but this region turns out 
to lead to a too large Z-exchange direct detection contribution. 
As a result this scenario is nothing but excluded. 

Next we looked at the possibility that the inert scalar asym¬ 
metry produced could still be at the origin of the observed 
DM relic density, which could be the case if it is reprocessed 
to a lighter specie, S, which satisfies the direct detection con¬ 
straints. We considered 2 scenarios where DM is made of 
a singlet odd under the Z 2 symmetry, a) Slow decay of the 
asymmetry into the (real or complex) singlet particle, occur¬ 
ring after S freeze-out. b) Reprocessing of the inert scalar 
doublet asymmetry into a S asymmetry through faster decays 
occurring before S freeze-out. Both possibilities can lead to 
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the observed relic density provided the interaction causing the 
decay is small enough to induce this decay at the right time. 

As most asymmetric DM scenarios, the framework we con¬ 
sider does not explain why the baryon and DM abundances are 
so similar. Our scenario trades this abundance coincidence for 
a coincidence between the mass of the proton, the mass of the 
inert states, the mass of the dark matter particle, and the val¬ 
ues of various couplings. Even if both abundances have same 
origin, these parameters must ’’cooperate” to lead to a DM 
abundance so close to the baryon one. Rather than providing 
a real explanation for the abundance coincidence, this scenario 
shows instead that the origin of the DM relic density could be 
of asymmetric origin, due to the generation of an inert scalar 
asymmetry related to the generation of a B — L asymmetry at 
high temperature. 
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In Eq. ( [ 15 ] ) and (29 1 , the effective cross section of the H 2 H 2 
SM SM coannihilations is given by 13 


yeq j 


+ j ■ 

(45) 

where g is the weak coupling constant. We neglected the X 5 
contribution, and the corrections due to the contributions pro¬ 
portional to (v^). 


Analytical resolution of the Boltzmann equations 


The Boltzmann equations given in Eqs. (|24|l-(26 1 do not in 
general have a simple analytical solution. However, in the 
case of very fast oscillations, like it is the case here, a good 
approximation consists in symmetrizing the equations for A^o 
and E^o, i.e. replacing Eqs. (251-(261 by 


ENRS-ERS, the IISN, the ERIA and an ULB-ARC and the 

BA^o 

-..■S'” ^ „ 

1 (Oov). 

Belgian Science Policy, lAP VTl 1. 

dz 

1 

t 

( 

J 

C 

2 zH 


BE^o 

2/^"* A 0 

1 (oov). 


dz 

zH ’l 

2 zH 


A^o 


i-rjO E .^0 


(46) 

(47) 


Appendix 

Rates and cross sections 


In Eq. (|^, the reaction density of the A ,5 scatterings for the 
r|+ (and similarly for ri**) is given by 




(39) 


where 


4 ^ 

'fcd = v^ f dxy/xKi{zy/x)a{ab ^ cd) . (40) 

647tz J 4 

with a{ab —i' cd) the reduced cross section. These are given 
by 


In this approximation, the solutions for A^o and E^o are of the 
form 

A^o(z) =/(z) cos[g(z)] , E^o(z) = //(z) sin[g(z)] . (48) 
Eurthermore, since we are interested in oscillations happening 


after the freeze-out, we can neglect in Eq. (|24|l. With these 


approximations, integrating from zew to z withthe initial con¬ 
ditions Aqg(z£w) = St|o(z£w) and Z{zew) = 0, the analytical 
solutions of the Boltzmann equations Eqs. are given 

by Eq. ( |48l l and 


Et |0 


(z) = y^A^oW-S2o(z) = /(z) , (49) 


with 






3X2 



fi^) = 

St|o (^Ew) 

(50) 


(41) 

1 + 5%“ (a 

27t ' 

Y-iy. 

y xj 

(42) 

g(z) = 

5m / z^ \ 

mz) Ulw J ■ 

(51) 


In the non-relativistic limit, the corresponding rate is given by 

(43) 


= ^ = {q^v) 


where 


3X2 


3X2 


= 32^ ’ 


H2 


16 Km 


H 2 


The abundance E^o decreases therefore monotonically until it 
reaches an asymptotical value given by 


E^o(z » Zew) = 


Srio i^Ew) 


^ + 2 H{z) z£^2-iio(Z£W) 


(52) 


Note that despite appearance, the denominator doesn’t depend 
on z, since sz/{Hzew) = ^'2y^MpiTEw/5K^- 
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